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Marsh Fritillary (Euphydryas aurinia) in the Czech Republic: 
monitoring, metapopulation structure, and conservation 
of an endangered butterfly 
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1. Introduction 


Hula, V., Konvicka, M., Pavlicko, A. & Zdenek, F. 2004: Marsh Fritillary (Euphy- 
dryas aurinia) in the Czech Republic: monitoring, metapopulation structure, and 
conservation of an endangered butterfly. — Entomol. Fennica 15: 231-241. 


Thirty colonies of the Marsh Fritillary (Euphydryas aurinia) have recently been 
recorded in Western Bohemia, Czech Republic. The majority of colonies are 
small, their total area are 1.5 km’. Small size and intensive grazing/mowing were 
positively associated with observed declines/extinctions, while abandonment 
threatens the colonies in the longer term. Short distances to nearest colonies buf- 
fered against declines. High colony turnover, asynchronous local dynamics 
pointed and the species’ biotope requirements all point to a dynamic metapopula- 
tion structure; patterns of connectivity revealed that there are several metapopu- 
lations within the region. Because conserving the species within its extant sites 
seems unsustainable in the long term, restoration of its habitats is proposed. 
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donment of marginal lands (Van Swaay & Warren 
1999). It is included in Annex II of the EEC/EU 


The Marsh Fritillary Euphydryas aurinia (Rot- 
temburg, 1775) (Nymphalidae: Melitaeinae) 
ranks among the most threatened butterflies in 
Europe. This inhabitant of unimproved marshy 
meadows has declined continentally due to land 
drainage, agricultural intensification, and aban- 


Habitat Directive, and legally protected in many 
European countries. The butterfly is highly selec- 
tive for both quality and spatial distribution of 
habitats (Porter 1984, Munguira et al. 1997, Le- 
wis & Hurford 1997). It persists in landscapes via 
metapopulation dynamics, forming spatially re- 


10 Kilometers 


stricted colonies interconnected by adult dis- 
persal (Warren et al. 1994, Wahlberg et al. 
2002a). The practices maintaining its habitats dif- 
fer regionally: British populations inhabit pas- 
tures (Asher et al. 2001), hay meadows are colo- 
nised in Central Europe (Anthes et al. 2003a, b, 
Konvicka et al. 2003), and woodland clearings 
are inhabited in Scandinavia (Wahlberg et al. 
2002b). In addition, the species varies geographi- 
cally in host plant use (e.g., Singer et al. 2002, 
Wahlberg et al. 2001, Wang et al. 2004), life his- 
tory (Eliasson & Shaw 2003) and genetic struc- 
ture of populations (Descimon et al. 2001). 

Until recently, little has been known about the 
ecology of E. aurinia in Central Europe. In the 
Czech Republic, a detailed survey of the species’ 
distribution has been carried out since 2000. It has 
revealed that the butterfly retreated from over 
60% of its 20th century range and became re- 
stricted to the region of Western Bohemia (Hula 
etal. 2002; see also Fig. 1). A study of larval habi- 
tat selection revealed a preference for large and 
clumped host plants growing at dry, nitrogen- 
poor and short-sward patches within humid 
meadows (Konvicka et al. 2003), and similar pat- 
terns apply to prealpine Bavaria (Anthes et al. 
2003a, b). 

Here, we focus on landscape-level patterns in 
distribution of Euphydryas aurinia in the Czech 
Republic. We describe the spatial distribution of 
colonies; estimate their size, report on manage- 
ment and currently faced threats; and analyse 
changes in colony sizes observed during a period 
of three years. We document that the species ex- 
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Fig. 1. Western Bohe- 
mia, with positions of ex- 
tant (and recently ex- 
tinct) colonies of 
Euphydryas aurinia. The 
symbols distinguish 
large colonies (>50 lar- 
val nests; empty circles), 
small to intermediate 
colonies (circles with 
dark centres), and colo- 
nies that had gone ex- 
tinct in 2001-2003 
(crossed circles). 


hibits a metapopulation structure, argue that clus- 
ters of adjoining colonies should be conserved as 
individual units, and that landscape-level restora- 
tion of habitats is necessary to safeguard the but- 
terfly in the long-term. 


2. Material and methods 
2.1. The butterfly 


Adult flight of E. aurinia lasts from the end of 
May until late June in Western Bohemia. Females 
lay clusters of eggs on leaves of their sole region- 
ally used host plant, Succisa pratensis (Moench, 
1794). Larvae hatch in July and feed communally 
in silken webs. They reach fourth instar in Sep- 
tember and weave more compact hibernation 
webs hidden near ground within S. pratensis ro- 
settes. They resume feeding shortly after 
snowmelt and feed first in loose groups and then 
solitarily on fresh S. pratensis leaves until pupa- 
tion. 


2.2. Survey of colonies 


Prior to 2000, only four colonies (numbers 7, 18, 
19, 20 in Fig. 1) of E. aurinia were known to local 
lepidopterists (Konvicka et al. 2003). We 
searched for additional colonies by visiting his- 
torical sites and checking other humid grasslands 
throughout the area of distribution of the butter- 
fly, using recollections of local naturalists, topo- 
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Table 1. Overview of data collected during monitoring of larval webs, and mark-recapture of adults, of 
Euphydryas aurinia in the Czech Republic, 2001-2003. Numbers of sites correspond with the numbers in Fig. 1. 
Larval abundances are expressed on a semiquantitative scale: 0; 1 web; 2: <5 webs; 3: <10; 4: <20; 5: <30; 6: 
<50; and 7: <100. ex = extinct, i.e. the site was visited, but no larval webs were found, n.a. = data not available. 











2001 2002 2003 
Site larvae larvae days marked recaptures larvae days marked recaptures 
1 4 4 25 559 400 7 5 298 156 
2 n.a. 4 n.a. — — 5 4 49 20 
3 3 3 23 229 105 3 20 254 228 
4 n.a. 4 23 228 202 4 21 369 497 
5 4 3 23 96 127 ex 20 136 189 
6 n.a. 3 23 86 40 3 20 138 178 
7 4 3 23 96 137 2 20 65 93 
8 n.a. n.a. n.a. — — 4 1 7 1 
9 2 ex 1 1 0 ex n.a. — — 
10 4 ex 1 7 3 ex 1 4 1 
11 n.a. 3 3 57 14 3 1 102 29 
12 3 ex 3 0 0 ex 1 0 0 
13 2 1 2 5 2 1 1 1 1 
14 5 4 1 43 13 5 3 103 22 
15 3 2 n.a — — 2 3 15 16 
16 1 ex n.a — — ex n.a — — 
17 4 3 2 40 10 ex 1 0 0 
18 4 1 n.a. — — 2 n.a. — — 
19 4 2 n.a. — — 2 n.a. — — 
20 3 ex n.a. — — ex n.a. — — 
21 n.a. 7 n.a. — — 7 n.a. — — 
22 n.a. 7 n.a. — — T n.a. — — 
23 n.a. n.a. n.a. — — 4 n.a. — — 
24 n.a. 6 n.a. — — 6 n.a. — — 
25 n.a. T n.a. — — 7 1 101 T 
26 n.a. 4 n.a. — — 4 1 46 10 
27 n.a. 2 n.a. — — 2 1 4 1 
28 n.a. n.a. n.a. — — 2 1 6 2 
29 n.a. 2 n.a. — — 2 1 8 8 
30 3 3 1 23 5 4 n.a. — — 





graphic maps, aerial photographs, and informa- 
tion of host plant occurrence from vegetation sur- 
veys. We thus visited some 100 sites and discov- 
ered 26 hitherto unknown colonies. The surveys 
took place mainly during late August/early Sep- 
tember, when the larval webs are most conspicu- 
ous, but some sites were discovered during adult 
period as well. 


2.3. Population size estimates 


During autumn visits, we zigzagged each site and 
recorded the abundance of larval webs on the fol- 


lowing ordinal scale: 0; 1; 2: <5 webs; 3: <10; 
4: <20; 5: <30; 6: <50; and 7: <100. 

To compare larval counts with adult numbers, 
we needed a method for estimating adult abun- 
dance across many sites. Thomas (1983) pro- 
posed performing quick estimates across multiple 
colonies, fitting them to estimates from a well- 
studied reference colony, and correcting the result 
for phenological effects of butterfly emergence. 
We selected one area with a high density of local 
colonies (hereafter “main area”) for detailed 
mark—release-recapture (MRR) survey. There 
were five colonies within 2 km’ (Fig. 1:3, 4, 5, 6, 
7), and an additional three colonies 3 km apart 
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In N (Jolly) 











In N (Craig) 
Fig. 2. Regression of Craig's vs. Jolly’s daily popula- 
tion size estimates obtained from main area sites. Re- 
gression based on 2002 and 2003 results; full line = 
fitted model, broken line: y = x. 


(Fig. 1: 1, 2, 8). The MRR was carried out in 2002 
(all colonies) and 2003 (the five colonies within 
2 km’). Marking proceeded on an every day basis, 
weather permitting, during the entire adult flight. 
We spent approximately equal time within each 
colony, marked all netted butterflies with unique 
numbers, recorded their sex, position of capture 
and wing wear (scale 1—4), and released them at 
capture points. In addition, we visited as many as 
possible of the more distant (“peripheral”) colo- 
nies, and spent there two to four person-hours per 
visit. Table 1 summarises the marking effort. 

We calculated two estimates of adult numbers 
for the main area sites. The first was the Jolly- 
Seber-Cormack [JSC] method (e.g. Jolly 1965), 
which gives one-day estimates for intensively 
sampled populations with births, deaths and in- 
terchange of individuals (Williams et al. 2001), 
using model A in the program Jolly (Pollock et al. 
1990). 

The second was the single-marking method of 
Craig (1953), which estimates population size 
from numbers of animals and numbers of capture 
events, assuming Poisson’s distribution of indi- 
viduals captured zero times, once, two times, etc. 
(Pollard 1977). We used only the latter method 
for peripheral colonies. 

The Craig method frequently overestimates 
true population sizes (Gall 1984). To correct for 
the bias, we regressed, separately for each study 
year, the (In-transformed) daily estimates from 
main area sites obtained by the two methods (Fig. 
2). We then used the regressions to convert the 
Craig estimates from peripheral sites into respec- 
tive JSC estimates. 

It remained to convert the JSC estimates from 
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peripheral sites to a mutually comparable pheno- 
logical date. We used colony 4 from the main area 
(the densest data available) as a template, and 
converted the JSC estimates from peripheral sites 
into the JSC estimate expectable at peak flight 
(Fig. 3). 

To estimate adult numbers for the colonies 
where no marking was done, we regressed, sepa- 
rately for each year, the (In-transformed) peak- 
flight JSC estimates against larval web counts for 
the colonies from which both estimates were 
available, and used the regressions to predict peak 
adult numbers for remaining colonies. 


2.4. Parameters of sites 


Individual colonies were characterised by (1) 
Area, defined as distinct biotope borders such as 
margins of forest or improved meadow; (2) Host 
plant abundance (ordinal scale, 1: <10, 2: <50, 
3: <100, 4: <500, 5: <1,000, and 6: above 1,000 
rosettes); (3) Management (None, Intensive pas- 
ture, Random [= light, or erratic] pasture, Mow- 
ing once a year, Mowing twice a year, and Tram- 
pling by game animals), recorded as a percentage 
of patch area (rounded to tens of per cent); and (4) 
Current threats, i.e., overgrazing, intensive hay- 
making, successional overgrown, afforestation, 
nutrient spillovers from adjoining fields, and land 
drainage, again in percentages of patch area af- 
fected. 

To express the connectivity of sites, we used 
two measures. The neighbourhood distance 
[“Neighbourhood”’] was simple centre-to-centre 
distance to the nearest colony. The landscape 
connectivity [“Connectivity”] took into account 
the distances of a given site to all other sites, 
weighting the distances by site areas. Following 
Moilanen & Nieminen (2002), 


S=) ae (1) 


where A, are areas of surrounding sites and dare 
centre-to-centre Euclidean distances from site j to 
k. The parameters a and ¢,_ scale the effect of dis- 
tance on migration and the effect of patch area on 
immigration, respectively. We estimated them 
from MRR data on dispersal, 2002, sexes com- 
bined, using the Virtual Migration model by 
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Fig. 3. Method used for recalculation of estimates of adult numbers from peripheral 
colonies to maximum potential adult numbers, year 2002, males. — a. Capture sex ra- 
tio fitted against days (straight line) indicate daily phenological states of the popula- 
tion. — b. Wing wear fitted against days (straight line) indicate daily phenological 
states of the population. In (a) and (b), black squares = Jolly’s one-day estimates from 
the well-studied colony 4 that were fitted using 3rd period moving mean; by setting 
maximum daily estimates as 100%, the resulting line gives a proportion of the maxi- 
mum population size present in particular days. The recalculation of estimates (for a 
and b) from peripheral colonies proceeds as follows: (i) Craig’s estimates are trans- 
ferred to JSC estimates, using Equation 5, above; (ii) capture sex ratio and wing wear 
recorded from a peripheral colony give proportions of maximum population size pres- 
ent there at marking; (iii) if the two proportions differ, their middle value is used; (iv) 
the estimated population size is then transferred to maximum population size, i.e. 


100%. 


Hanski et al. (2000). The values were a = 1.5 and 
Gm = 0.11. Details of the Virtual Migration ap- 
proach will be presented elsewhere. 


2.5. Predictors of decrease and extinctions 


We constructed multiple regression models of de- 
crease/extinction against site parameters, using 
generalised linear models with binary dependent 
variables (decrease vs. stability and/or increase, 
and extinction vs. persistence) and logit-link 
function. For “decreasing”, we considered the 
colonies where the abundance of larval nests 
dropped by 2 points at the ordinal scale used for 


nests counts (see section 2.3.) during the three 
years. Percentual predictors (i.e., those referring 
to management) were arcsine transformed prior 
the analyses. Following single-term regressions 
for all predictors, we used stepwise forward se- 
lection to select the best-fitting model, comparing 
model fits against respective null models using 
the Akaike information criterion (AIC) and F- 
tests. We used S-plus 2000 (S-plus 1999) for the 
computations. 


3. Results 


Fig. 1 shows the locations of 30 extant (and re- 
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cently extinct) colonies of E. aurinia in Western 
Bohemia. The majority of the colonies belong to 
distinct clusters, but some colonies are separated 
from the nearest colony by >10 km of unoccupied 
areas. 


3.1. Estimates of population size 


We marked 1595 adults in 2002 and 1793 in 
2003, and obtained 1097 and 1558 recaptures 
(Table 1). For the main area (sites 1, 3, 4, 5, 6 and 
7 in 2002 and 3, 4, 5, 6 and 7 in 2003), regressing 
of the daily Craig and Jolly estimates gave the re- 
lationships: 











In N sc = 0.72 (0.067 S. E.) * In (N, oa) 

+ 0.96 (£0.340 S. E.) (year 2002) (2) 
and 

In N sc = 0.86 (+0.104 S. E.) * In (N, oa) 

— 0.03 (+0.498 S. E.) (year 2003) (3) 





Both regressions were significant (2002: F, 12 = 
113.3, p <0.001, R = 0.49; 2003: Fia” 680, 
p <0.001, R? = 0.61) and neither their slopes nor 
their intercepts differed (slope: t, = —1.75, two- 
tailed p = 0.08; intercept: t, = —1.84, two-tailed 
p=0.07). We hence constructed the common re- 
gression, 

In N sc = 0.76 0.057 S. E.) * In Wenig) 

+ 0.64 (+0.284 S. E.) (4) 








Plotting it against y = x line (Fig. 2) showed that 
Craig’s method underestimated population size 
for low adult numbers (<100 individuals), and 
overestimated it for large populations (>1,000 in- 
dividuals). 

Fig. 3 illustrates the method used for re-calcu- 
lating one-day adult estimates to numbers in peak 
flight. Use of the method was justified, because 
the male to female ratio decreased with the season 
(logit regression, 2002: %? = 38.3, df = 1, 
p <0.0001; 2003: x? = 26.4, d.f. = 1, p <0. oe 
and wing wear gradually increased (2002: F, „= 
344.3, p <0.001; 2003: F, ,, = 90.3, p <0.001). 

The regressions of the O SC-estimated) peak 
adult numbers against larval counts were 
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In (Nigeoma,) = 0-78 (£0.29 S. E.) * (nest count) 
+ 2.69 (+ 0.873 S. E.) (year 2002) (5) 





and 





In (Ncn) = 0-58 (+0.100 S. E.) * (nest count) 
+2.58 (+ 0.467 S. E.) (year 2003) (6) 





They were significant (2002: F, , = 7.2 »P <0.05, 
R = 0.44, 2003: Hg aoFs p<0.001, R =0.69) 
but differed in slopes (slopes: £, = 2.39, p <0.05, 
intercepts: t,, = 1.89, p = 0.07). Hence, we could 
not use a single relationship for recalculating lar- 
val and adult counts. 

Using equations 6 and 7, we estimated the to- 
tal Czech population of E. aurinia. For 2003, i.e. 
the year with the most complete survey, the esti- 
mate was ~ 5,000 using the 2003 parameters and 
=~ 8,000 using mid-values between parameters 
from 2002 and 2003 (slope: 0.65, intercept: 2.62). 
Although such estimates are crude and total 
emergence exceeds the numbers in peak flight, 
they illustrate how small the Czech population of 
E. aurinia in fact is. Regarding year-to-year vari- 
ation, the estimates from 11 sites surveyed each 
year gave annual estimates ~ 2900 (2001), ~ 1700 
(2002) and ~ 1500 (2003), with the coefficient of 
variation 0.28. 


3.2. Parameters of colonies 


Out of the 30 extant colonies, 23 contained <10 
larval nests (i.e., ca. 70 adults in 2002 and ca. 150 
adults in 2003) whereas four contained >100 
nests or >1,000 adults. Areas of 23 colonies were 
smaller than five hectares, whereas two colonies 
exceeded 25 hectares. The summed inhabited 
area was 1.44 kin’. 

If management of the colonies is expressed 
proportionally to their total area, than two thirds 
of the area is neglected land (neither mown nor 
grazed, while 30% of the total area is either 
grazed or mowed (Fig. 4a). Notably, although 
trampling of game animals affect only small pro- 
portions the sites, the host plants of E. aurinia 
tend to concentrate at these disturbed patches (cf. 
Konvicka et al. 2003). 

Twelve of the colonies (40% of total) are 
threatened by succession or by afforestation. On 
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Fig. 4. Overview of recent (2003) management, and 
threats, of the 30 sites inhabited in 2000-2003 by 
Euphydryas aurinia in Western Bohemia, including 
sites where the butterfly had gone extinct during the 
period. — a. Management, recalculated to percentage 
cover of individual sites. — b. Threats, related to the 
percentage of thus affected individual sites. 


the contrary, ten sites are threatened by over- 
management: they are mown or grazed too inten- 
sively, affected by past drainage, or by leakage of 
nutrients from adjoining farms (Fig. 4b). Related 
to total inhabited area, the picture is slightly more 
optimistic, as one-third of the total is not threat- 
ened recently, while one-third is threatened by 
succession/afforestation and one-third by too-in- 
tensive farming. 


3.3. Evidence for metapopulation dynamics 


Two prerequisites of metapopulation dynamics 
are asynchronous within-site dynamics and lim- 
ited, but existing, dispersal (Hanski 1999). Re- 
garding the former, three out of 24 colonies sur- 
veyed for >1 year increased in numbers, twelve 
decreased and seven went extinct (Fig. 1). Re- 
garding dispersal, it suffices to state that out of 
600 butterflies recaptured at least once in 2002, 
378 did not leave their sites of original capture. 
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The mean movement was 151 (SD = 537.8) me- 
ters, 22 butterflies crossed distances >1 km 
(males: 18, females: 4), and one female was cap- 
tured after dispersing from site 14 to 11 (8.5 km). 
Hence, some individuals were capable of consid- 
erable inter-colony movements. 


3.4. Correlates of decrease and extinction 


Decreases in colony size (Table 2) were less 
likely at large sites that contained high numbers 
of rosettes of S. pratensis and were mown once a 
year, and more likely at intensively grazed sites. 
They were more likely the longer the distance to 
the nearest colony, but did not depend on land- 
scape connectivity (S). The multiple-regression 
model accounted for 51.7% of the variation and 
contained hostplant abundance (positively corre- 
lated with area, Spearman’s r = 0.71, p <0.001) 
and the distance to the nearest colony. 
Extinctions (Table 3) were less likely in colo- 
nies with high abundance of host plants that were 
mown once a year. Positive predictors of extinc- 
tion included grazing (both intensive and erratic) 
and intensive (twice-a-year) mowing. Host plant 
abundance, which was the only predictor of the 
final model, explained 47.2% of the variation. 


4. Discussion 


Although increased survey intensity consider- 
ably increased the number of known sites of £. 
aurinia in the Czech Republic, most of the colo- 
nies are small, the entire flight area does not ex- 
ceed 1.5km’, and the butterfly remains seriously 
threatened. 


4.1. Monitoring 


Our estimates of the population size from periph- 
eral localities relied on comparisons with inten- 
sively sampled sites. The fact that the relation- 
ships between estimates obtained by the quick 
Craig’s and relatively laborious JSC methods 
were consistent across two successive years may 
greatly aid to future monitoring. Naturally, our 
approach relied on several simplifying assump- 
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Table 2. Factors associated with decreases of colony size (expressed as a decline of larval web abundance by 
>2 ordinal points used for larval surveys). GLM regressions with binarily coded dependent variable (decrease vs. 
stability of increase). “Direction” indicates positive or negative relationship, F-tests refer to comparisons with null 
model. Only colonies followed for at least two years were analysed. * p <0.05; ** p <0.01; *** p <0.001. 








Direction Resid. df Deviance Cp F p 
deviance 

Null model 35.6 26 37.4 
Single-term models 
Area — 28.6 1, 24 7.0 34.1 3.8 * 
Connectivity 33.5 1,24 2.1 40.0 0.4 n.s. 
Neighbourhood + 30.4 1,24 5.2 34.5 5.1 7 
No management 35.6 1,24 0.0 39.9 0.0 n.s. 
Mowing once — 21:7 1,24 7.9 30.9 9.8 sd 
Mowing twice 35.5 1,24 0.1 39.8 0.1 n.s. 
Intensive pasture + 31.3 1, 24 4.3 35.3 4.3 * 
Random pasture 32.8 1,24 2.8 37.0 2.8 n.s. 
Trampling 34.7 1,24 0.9 39.1 0.8 n.s. 
Host plant abundance — 20.9 1, 24 14.7 24.2 IT-7 aa 
Best model 
Hostplant + neighbourhood 17.2 2, 23 18.4 22.1 11.3 a 





Table 3. Factors predicting extinctions of colonies. GLM regressions with binarily coded dependent variable (ex- 
tinct vs. not extinct). F-tests refer to comparisons with null model, only colonies followed for at least two years 
were used in the analysis. “Hostplant’ was the only variable that entered the final model. * p <0.05; ** p <0.01; 








*** o <0.001. 

Direction Resid. df Deviance Cp F p 

deviance 

Null model 30.9 26 32.8 
Single-term models 
Area — 25.8 2, 24 5.1 31.4 2.7 + 
Connectivity 27.7 2,24 3.2 34.1 1.5 n.s. 
Neighbourhood 30.9 1,25 0.1 35.2 0.1 n.s. 
No management 29.2 1,25 1.8 33.6 1.6 n.s. 
Mowing once - 25.9 1,25 5.0 29.1 6.3 z 
Mowing twice 30.8 1,25 0.1 35.2 0.1 n.s. 
Intensive pasture + 25.0 1,25 5.9 29.0 5.9 gj 
Random pasture + 26.0 1,25 4.9 30.2 4.7 i 
Trampling 30.6 1,25 0:3 34.9 0.3 n.s. 
Hostplant abundance — 16.3 1,25 14.6 20.7 13.4 = 





tions, including reliability of the JSC method, 
monotonous increase of wing wear and monoto- 
nous decrease of sex ratio (i.e., protandry) during 
season. However, the JSC method can always be 
substituted by more sophisticated abundance esti- 
mates (cf. Schtickzelle et al. 2003), and both 
protandry and increasing wing wear are com- 
monplace in butterflies with discrete generations 


(e.g., Fagerstrom & Wiklund 1982). 

The regressions of adult numbers against lar- 
val counts differed between the years. This might 
be due to low sample size in 2002, when larval 
and adult counts were available for just eleven 
colonies. In addition, changing management of 
the sites could distort any relationships between 
adult and larval numbers. We observed this at site 
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5, in 2003, where females were quite abundant 
(captures: 61, NJSC_, = 92) but all nests were 
destroyed by mowing in August. Hence, caution 
is required for using larval abundances to predict 
adult numbers, and vice versa. Regardless, the 
counts of larval webs sufficed to detect relative 
population changes, and thus remain efficient for 
cheap and quick monitoring. 


4.2. Metapopulation structure 


The asynchronous dynamics of local colonies 
and limited dispersal pointed to a metapopulation 
structure similar to that found in Bavaria (Anthes 
et al. 2003a, b) and other regions (Warren 1994, 
Wahlberg et al. 2002b, Wang et al. 2003, 2004). 
This similarity among distant regions supports 
the notion that metapopulation dynamics is an in- 
trinsic trait of this species (Wahlberg 2000). In 
agreement with metapopulation theory, small 
area and isolation of patches adversely affected 
colony survival. However, we did not find any ef- 
fect of overall landscape connectivity. This indi- 
cates that dynamics of individual colonies were 
influenced by the situation in their close vicinity 
rather than by the situation in the entire area, and 
that there are several mutually independent 
metapopulation systems. Since the proximity to 
other colonies buffered against decline and ex- 
tinction, these systems should be managed as in- 
dividual conservation units. Joyce & Pullin 
(2003) observed a similar situation in Britain, 
where no genetic differentiation existed within 
clusters of E. aurinia colonies, but isolation by 
distance was evident among the clusters. 

The area surveyed in this study was relatively 
sparsely occupied compared to, e.g., the meta- 
population of Melitaea cinxia on the Aland Is- 
lands, Finland (Hanski 1999). Although there 
still is a possibility that we failed to discover some 
colonies, it is unlikely that we missed many of 
them, or that some of the undetected colonies are 
large. The entire area suffered massive land 
drainage and meadow improvement schemes in 
the 1970s and 1980s, and the remaining humid 
meadows form tiny, and easily distinguished, is- 
lets in an otherwise unsuitable landscape. 

Contrary to most metapopulation studies, we 
did not consider the existence of empty patches. 
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Dennis et al. (2003) argued that while it is easy to 
spot a patch if the organism is present, without it 
one can never be certain whether a site contains 
all necessary resources. For E. aurinia, the simple 
presence of the host plant does not suffice to de- 
fine habitat, as its larvae prefer a specific vegeta- 
tion structure (Konvicka et al. 2003) and adults 
require rich nectar supply and sheltered spots uti- 
lised for perching (unpublished observations). 
Separation of mating and egg-laying patches was 
reported for at least one checkerspot species 
(Schwarzwalder et al. 1997). The uncertainty re- 
garding empty patches explains why we did not 
observe any (re)colonisation. To observe it, we 
should have observed extinction first, and ob- 
serving both events at one spot was very unlikely 
during our short-term study. 


4.3. Site management 
and conservation guidelines 


Anthes et al. (2003a, b) argued that joint effects of 
quality and spatial distribution of habitats affects 
the persistence of E. aurinia in Bavaria. Our ob- 
servation that intensive grazing and mowing had 
deleterious effects on colony persistence corrob- 
orates this claim. Intensive grazing reduces nectar 
availability (Dolek & Geyer 1997, Anthes et al. 
2003a) and host plant size (Konvicka et al. 2003), 
whereas late-season mowing destroys larval 
webs. Less expectable were the deleterious ef- 
fects of light grazing (Table 3). We observed the 
effect at two fenced streamsides adjoining im- 
proved pastures (sites 10 and 17), where only a 
few days of cattle presence sufficed to destroy 
host plants with egg batches. However, both colo- 
nies were very small prior to this event. 

On the other hand, prolonged abandonment 
leads to encroachment of sites by tall vegetation, 
shrubs and trees. These contradicting effects 
point to the paradox that conditions beneficial for 
the butterfly develop shortly after cessation of 
mowing or grazing, but persist for a relatively 
short time. In the past, traditional farming had 
maintained fine-grained landscape mosaics of 
varying managements, which allowed the butter- 
fly to persist by tracking transitorily neglected 
patches (cf. Wahlberg et al. 2002b). However, 
such landscape mosaics no longer exist. Al- 
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though we might try to mimic them within indi- 
vidual colonies (Konvicka et al. 2003), the high 
colony turnover indicates that restricting conser- 
vation efforts to individual colonies would be un- 
sustainable in the long term. 

Extending the recently inhabited area offers 
the only long-term choice. The methods available 
include re-establishment of mowing at aban- 
doned sites, and diminishing grazing/mowing in- 
tensity at improved meadows and pastures 
(Barnett & Warren 1995). At present, both the 
space and resources for a landscape-scale restora- 
tion are available. The butterfly inhabits regions 
unsuitable for intensive agriculture, where farm- 
ing already depends on government subsidies. 
The total sums spent annually are difficult to as- 
certain, but it is illustrative that the single farmer 
who manages our main area receives annual sub- 
sidies for 300 hectares, i.e., for more land than is 
inhabited by the species at present. Therefore, 
managing the area for the butterfly should not in- 
cur any additional cost to taxpayers. 

A species action plan for E. aurinia, now un- 
der preparation in the Czech Republic, advocates 
the following modifications of the conditions 
qualifying local farmers for subsidies: execution 
of hay harvests in a patchy manner with portions 
of meadows left uncut each year, lower stocking 
of pastures, and not fertilising meadows. It also 
proposes cessation of governmental support for 
afforestation, since wooded conditions restrict 
adult movements (Norberg et al. 2002). The ulti- 
mate goal is to enhance the inhabited area five- 
fold, which seems realistic given the extent of 
subsidised (or abandoned) grasslands. 
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